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Abstract 

The paper gives the polarization of the tau lepton in the semileptonic B 
decays with respect to the direction of the virtual W boson. The result is 
given including the nonperturbative HQET corrections. The perturbative QCD 
corrections are probably negligible as suggested by the existing results for the 
longitudinal polarization of the charged lepton (Jezabek and Urban, 1998). 
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1 Introduction 



The interest in semileptonic B decays is currently increasing as the B factory in 
KEK is scheduled to begin to collect data later this year. This domain of physics 
is likely to upgrade our knowledge on the Standard Model parameters as well as to 
provide tests on its validity. The semileptonic B decays can contribute to the former 
as their theoretical description is now far more successful than that of the hadronic 
processes |]]- |J. 

The polarization of the charged lepton does not depend on the Cabibbo-Kobayashi- 
Maskawa matrix element and so it can be instrumental in finding the quark masses. 
The longitudinal polarization of the tauon including first order perturbative QCD 
corrections has been found analytically || by taking the analytical decay width for 
the unpolarized case |7j and calculating the width for a negative polarization. Then 
the result can be integrated to give polarized tau energy distribution. The method 
used in that calculation can easily be modified to give other polarizations. This fact 
matters insomuch that experimentally it is the polarization along the intermediating 
W boson direction that is easier to measure ||, see also 0. The reason is that the 
direction of r lepton can be determined at B factories with rather poor accuracy. On 
the other hand the direction of W is opposite to the direction of hadrons in semilep- 
tonic B decays. The latter can be well measured at least for the exclusive B — > Drv T 
and B — > D*tu t channels which probably contribute the dominant contribution to 
the inclusive decay rate. 

The present calculation includes tree level and HQET corrections only. We are 
unable to calculate perturbative QCD corrections because the analytic structure of 
expressions is far more complicated than in the case of the longitudinal polarization 0. 
However, indications exist that the effects of perturbative QCD corrections on r po- 
larization are negligible. In particular, the above-mentioned longitudinal polarization 
does not change visibly after the first-order perturbative corrections have been in- 
cluded either in the rest frame of the W boson [|10J or that of the decaying quark 

The paper is broken up into four sections. In Sec. 2, kinematical variables are 
introduced. Sees. 3 and 4 are to explain the method used in the calculation and then 
the results are shown is Sec. 5. In Appendix A some details of HQET calculations are 
explained including the discussion of singularity problems. Such problems were also 



encountered in [11] where a method was proposed to eliminate them. 



2 Kinematical variables 

In this section we define the kinematical variables used throughout the article as well 
as their boundaries. The calculation is performed in the rest frame of the decaying B 
meson, which coincides with that of the b quark at the tree level in the parton model. 
The four-momenta of the particles are denoted as following: Q for the b quark, q 
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for the c quark, W for the virtual W boson, r for the charged lepton, and v for the 
corresponding antineutrino. All the particles are assumed to be on-shell so that their 
squared four-momenta equal their masses: 

Q 2 =m 2 , q 2 = m 2 c , r 2 = m 2 T , v 2 = . (1) 

The employed variables are scaled to the units of the decaying quark mass m&: 

m\ m 2 T 2E T W 2 2E V 

P = — j , V = — , y = , t = — , x = . (2) 

mi mi mi, mi m^ 

Henceforth we scale all quantities so that m 2 = Q 2 = 1. The charged lepton is 
described by the light-cone variables: 

T ± = l -(y±^y 2 -Ar ] y (3) 

The W boson is characterized likewise: 

w = ^(l+t-p), (4) 



W3 = yWo-t, (5) 
w± = w ± w 3 . (6) 

The phase space is defined by the ranges of the kinematical variables: 

2y/rj <y<l + v-p = ym, (7) 

tmin = T_Yl - Y~rj ~ * - T+ ( X ~~ 1 ^ = tmax ■ 

The limits above are obtained within the parton model approximation. They change 
if we allow for Fermi motion, which we must in order to be able to discuss the HQET 
corrections to the decay widths |TI]-[]I5]. Also, contrary to the parton model case, 



the energy of neutrino can vary within limits which depend in a non-trivial manner 
on the values of the variables y, t. The details of the subject were discussed in |T2 



so we will only state here that the integrations involving delta functions and their 
derivatives have the effect of confining the range of the variables y, t to that of the 
parton model. 

3 Polarization evaluation 

The polarization is found by evaluating the unpolarized decay width and any of the 
two corresponding to a definite polarization, according to the definition, 

r+ - r~ r~ 

P=- _ = l-2— , (9) 

r+ + r- r ' y ! 
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where r = T + + r~. The calculation of the polarized width is structured after the 
manner of that which has yielded the longitudinal polarization ||. Thus in the rest 
frame of the decaying quark, one can decompose: 

s = AQ + BW . (10) 

The coefficients A, B can be evaluated using the conditions defining the polarization 
four-vector s, which reduce to the following when the parton model value of the 
neutrino energy is assumed: 

*Z = T , t + V , (11) 

't(y-y-)(y+-y) 



t$ = ± = , (12) 

y/t(v-v-)\y+-v) 

where the superscripts at A, B denote the polarization of the lepton,while 

y ± = (l + ri/t)w± . (13) 

This observation is made relevant by the fact that the decay width for a definite 
polarization of the charged lepton is gotten from the analogous expression for the 
unpolarized case, 

dT = G 2 F M^\V CKM \ 2 M^dn 3 (Q; q, r, z/)/vr 5 (14) 

where 

M u %(r) = q- tQ ■ v , (15) 
by formally replacing the lepton four-momentum by the following four-vector K: 

K = t — m T s . (16) 

Then we obtain, 

Mli= l -MZ{K = T-ms)= l -{q-K){Q.v) . (17) 

Although the expressions above are written for the Born approximation, corrections 
received by the hadronic tensor obviously do not alter this scheme, so that we can 
apply it to the HQET calculations, too. However, the coefficients in the decomposition 
( |i~0l) need to be re-evaluated, taking into account the Fermi motion and working in the 
rest frame of the decaying meson. For a derivation of these cf. Appendix A. Applying 
now the representation (|i~0|) of the polarization four-vector s we readily obtain the 
following useful formula for the matrix element with the lepton polarized: 

M ± = l -M un {r) T , = ^ = -- x9 \yM un {W) - (t + V )M un (Q)] . 
2 2 ^y{t + r ] ){x + y)-yH~{t + 7 ] y 

(18) 
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The above expression is valid for the HQET corrections, too. The first term on the 
right hand side of can be calculated immediately once we know the result for 



the unpolarized case. Then the other terms require the formal replacement of the 
four- momenta r — > W and r — > Q in the argument. 



4 Evaluation of HQET corrections 

Using the operator expansion technique, one can obtain corrections to the decay 
widths of heavy hadrons which effectively lead to new terms in the hadronic tensor 
appearing in the triple differential decay width, 

dT WGl C^ . (19) 



dxdtdy 2tt 3 

The hadronic tensor W, related to an inclusive decay of a beautiful hadron H b , 

VIV = (27r) 3 5> 4 G^ ft -q-px) < H b (v,s)\tf\X >< X\J c u \H b (v, s) > (20) 
x 

can be expanded in the form 

= -9^W X + v^v v W 2 - ie^pV^Ws + q^q v W A + (g M u„ + v^Ws . (21) 

The form factors W n can be determined by using the relation between the tensor W 
and the matrix element of the transition operator 

= -iJ d 4 xe-^T[tf(x)r u (0)] , (22) 

which is 

W, lu = --Im< H b \T^\H b > . (23) 

The coefficients W n of (|21| ) have all been found elsewhere, see eg. |T2) for a complete 
list. Then the distribution (fl9|) can be schematically cast in the following form: 



<n f x S[x - x ) + f 2 S'(x - x ) + f 3 S"(x - x ) , (24) 



dxdtdy 
where 

x = 1 + t - p - y (25) 

is the value of the neutrino energy in the parton model kinematics. The triple dif- 
ferential distribution must be integrated over the neutrino energy to give meaningful 
results. The final lepton energy distribution obtained on two subsequent integrations 
may be trusted except for the endpoint region where the operator product expansion 
fails. In the present paper we give the double differential distribution so that the 
lepton energy distribution has to be obtained numerically. The calculation does not 
show any features unfamiliar from the cases of the other known polarizations. 
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5 Results 



5.1 Double differential distribution 

The polarized distribution can be written in the form, 

i dr± i 



r dy dt 



2 



± (F + F + 



where 



r 



G\m\ 

1927T 3 



CKM\ 



(26) 



(27) 



The first term on the right hand side of Eq. (|26|) stands for the unpolarized distribution, 
12H ) Eq.(30). Here we will only give the new other term: 2 



given e.g. m 
F = 

where 



fi 

h 
h 



W [6/i + K b W (/ 2 + f 3 W 2 + |/ 4 >V 4 ) + G b (/ 5 + f t W 



(28) 



-yt{l + p - 77) + 2yt 2 - y(l + p?/ - 2p + p 2 + 77) - y 2 t + y 2 (l - p) 

+ t(l + 2p7/-p 2 )+t 2 (2p-7/)-t 3 -p 2 77 + 77 , 

-8yt(l - p + 77) - 16yt 2 + 83/(1 + pr] - 2p + p 2 - r]) + QyH 

- 6y 2 {l - p) + 2t(l - 6pr/ - 4p + 3p 2 + 12r/) + 2t 2 (8 - 6p + 3r/) 
+ 6t 3 - 8p?7 + 6p 2 r] + 2rj + 8rf , 

yt(-lOpr) + Qpr/ 2 + Up 2 r] + 9pV - 6p 3 r? + 2rj + 9r/ 2 - Arf) 
+ yt 2 (l - 9pr/ 2 - 5p + 18p 2 r/ + 7p 2 - 3p 3 + 677 - 13r/ 2 ) 
+ 7//; 3 (l - 18pr/ - 2p + 9p 2 - 14t/ + 37/ 2 ) - yt\5 + 9p - 67/) 
+ 3?/t 5 + y(-5pf] 2 + Apr] 3 + 7pV - 3pV + r/ 2 + 4t/ 3 ) 
+ y 2 t(l + Upr] + 2pr] 2 - 3p - p 2 r/ + 3p 2 - p 3 - 13t/ + 8t/ 2 ) 
+ y 2 t 2 (-6 + 5pr7 + 6p + 15// - r/ 2 ) + y 2 t 3 {7 + 3p - 3r/) - 2?/ 2 t 4 
+ 7/ 2 r/(-3p + 8pr/ + 3p 2 - p 2 r] - p 3 - 7r] + 1) + ?/ 3 t(l - p) 2 

- y 3 * 3 + y 3 r]{l - p) 2 , 

y 2 t(-Qpr] 2 - 2pV + 6pV + 4p 3 r/ 3 - 3p 4 r/ 2 + 3r/ 2 - 2r/ 3 ) 



(29) 



(30) 



y 3 t 2 r) 



(31) 



+ y 2 t 2 (-6pr) - 2pr/ 3 - 6p 2 r/ 2 - 6p 2 r/ 3 + 6p 3 r/ + 12p 3 r/ 2 - 3p 4 r/ + 3r/ 



6r/ z 



+ y 2 t 3 (l - Qprf + 4pr/ 3 - 2p - 6p 2 r/ - 18pV + 12p 3 r/ + 2p 3 - p 4 - 677 + 2r/ 3 ) 

+ y 2 t 4 (-2 - Qpr] + 12pr/ 2 - 18p 2 r/ - 2p 2 + 4p 3 + 6r/ 2 - r/ 3 ) 

+ y 2 t 5 {12pr] - 2p - 6p 2 + 677 - 3r/ 2 ) + y 2 t 6 (2 + 4p - 3r/) - y Y 

+ y 2 (-2pr] 3 + 2p 3 rf - p A rf + rf) + y 3 t(8pr] + prf + 2pr/ 3 - 12p 2 r/ + pV 



+3p 2 r/ 3 + 8p a r] - p i rf - 2p 4 r/ - 2r) - rf + 3r/ d 



-,3^2 



2 A FORTRAN code for this formula is available from piotr@charm.phys.us.edu.pl 
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+ y 3 t 2 (-l - pr] + 6pr] 2 - 3pr] 3 + Ap - p 2 p + 9p 2 p 2 - 6p 2 + p 3 p + 4p 3 

-p 4 + 77 + 9p 2 - 3p 3 ) 
+ yH 3 (l + 6pr] - llprf - p + 9p 2 p - p 2 + p 3 + 9?7 - llr] 2 + rf) 
+ y 3 t 4 (3 - 13pp + 2p + 3p 2 - 13p + Arj 2 ) + y¥(-5 - 5p + 5p) + 2y¥ 
+ y 3 {Apif + pi] 3 - 6pV + p 2 i] 3 + Ap 3 if - p 3 if - p 4 rf - if - if) 
+ yH(-Qpi] + 2pi] 2 + Qp 2 i] + pV - 2p 3 p + 2i] - 3i] 2 ) 
+ y 4 t 2 (l + Apr] + pi] 2 - 3p + 2p 2 i] + 3p 2 - p 3 - 6p + 3p 2 ) 
+ y 4 t 3 (-3 + 2pp + 2p + p 2 + 6?7 - p 2 ) + y 4 t 4 (3 + p - 2p) 

- y 4 t 5 + y\-3pif + 3pV - p 3 if + if) , (32) 
/ 5 = Ayt{3 + 5p - 5?]) - 40yt 2 + 4y(l + 5pp - 6p + 5p 2 + 77) + 10y 2 t 

- 2y 2 (l - hp) - 6t(l + 10pr] - hp 2 ) - 2t 2 (A + 30p - 15p) + 30t 3 

+ 30p 2 p - 6i] + 8i] 2 , (33) 
/ 6 = yt(-18pp-2pp 2 + 6p 2 p-15pV + 10p 3 p + 2p + 17p 2 -4p 3 ) 
+ yt 2 {l - 16pi] + 15pp 2 - 9p - 30p 2 i] + 3p 2 + 5p 3 + IO77 - if) 
+ yt 3 (l + 30pr] - lOp - 15p 2 + IO77 - hrf) + yt 4 (7 + 15p - IO77) 

- 5yt 5 + y(-9pi] 2 + Apr] 3 + 3pV + 5p'V + p 2 + 8p 3 ) 

+ y 2 t(-l + 30pp - 10p?7 2 + 7p + 5p 2 p - lip 2 + 5p 3 - lip - 4p 2 ) 
+ y 2 t 2 (-2 - 25pp + 18p - 19p + 5p 2 ) - 15y 2 t 3 (l + p - 77) 
+ lOy 2 t 4 + y 2 (7pi] + \2prf - \\p 2 r] + 5p 2 p 2 + 5p 3 p - 77 - 9p 2 ) 
+ y 3 t(l - 6p + 5p 2 + 877) + y 3 t 2 {8 - 5p) - 5yH 3 

+ y 3 {-Qpr] + 5p 2 r] + 77) , (34) 



F± 
where 
hi,± 



h 



2,± 



l 3,± 



v/p W± { [tf 6 (h lt± + h^ ± W ± 2 ) + G b h 3 , ± ] 5(z±) + K b h 4 , ± 5'(z ± )} , (35) 



AyH 



12y 2 r] 



3y 2 t 



y 2 rj 



—8ytr] + 8yr] 

- 4<r± (6j/f - 2yr] 

- 8<r| (3j/t - 2/77 - 5y 2 + At + Arj) 
2a± y (8tp 2 + 16t 2 p + 8t 3 - Aytrf 



2y 3 t - 2y 3 r] - 2y 4 - Vbtt\ - 16p 2 



y 3 -8 



trj - At 2 - 4p 2 ) 



- 16a 3 ± (3y + t + r]) , 
Ayt 3 - 8y 2 tr] - 6y 2 t 2 
y 3 t 2 + y 3 rf + yH + y 4 r]) + Aa\ y {-Atrf - 8t 2 r] - At 3 



2^2 



2y z V 



(36) 



-12^77 - 8yt 2 - Ayr] 2 + 2y 2 tr] - y 2 t 2 + 3y 2 rf + 3y 3 t + 3y 3 r]) 

+ 8<r| y (-4*77 - 2t 2 - 2p 2 + 4ytp + yt 2 + 3yp 2 + 3y 2 t + 3y 2 r]) 

+ 16^ 1/ + t 2 + rf + yt + yrj) , (37) 

-16yt77 - 8yt + 8yt 2 + 8yr] + 8yp 2 - 8y 2 t + 8y 2 r] + 2At 2 - 24p 2 

+ 4cr± (-lOyt - 14yp - 5y 2 t + 5y 2 p - Ay 2 + 5y 3 + 12* - 4t 2 + 12p 
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h 



4,± 



+ At] 2 ) + 16a| (5yrj -2y + 5y 2 - 9t - 9r?) + 80a| (y + t + rj) 
4o± (%tr/ - %t 2 + 6y 2 t + 2y 2 r] + yH - y 3 7] - y A - 8^ - 8t 2 ) 
+ 8«4 (%t + %?7 + y 2 t - 3y 2 r? - 3y 3 + ttr) + At 2 ) 
- 16o-| (yt + 3yr/ + 3y 2 - 2t - 2r]) - 32o^ (y + t + rj) , 



and 



W4 



1 



y(t + r/)(2or± + y) - y 2 t - {t + V ) 2 
1 



^y(t + rj)(x Q +y)-yH-(t + rj) 2 
-- (t - ti)/(2t±) , z± = l + t-p-y-2a± 



(38) 

(39) 

(40) 

(41) 
(42) 



The parameters Kb, Gb, representing the kinetic energy and the chromomagnetic en- 
ergy, are defined according to |H 



5.2 Lepton energy distribution 

As regards the HQET correction terms, we only give the energy distribution in the 
form of a diagram evaluated numerically. Beneath we also give the Born level ap- 
proximation analytically The analytic formulae for the polarized distribution can be 
simplified if we split the kinematical range of y into two parts, separated by the value 
of the charged lepton energy where the virtual W boson can stay at rest. This value 
is 

y w = l-Jp+ T ^-_ . (43) 



In the formulae below, the superscripts A, B refer to the appropriate regimes: 



y <yw 
y>yw 



region A , 
region B . 



The energy distribution of polarized r lepton reads, 

dY ± n 



dy 



i2r n 



-f(y)±Af(y) 



(44) 
(45) 



(46) 



The function f(y) represents the unpolarized case, 

f(y) = ^C 2 vV-4r/ {C [y 2 - Sy(l + V ) + 877] + (3y - 677) (2 - y)} 

with 



c = i- 



p 

l + rj-y 



(47) 
(48) 
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The function Af(y) reads, 



Af(y) = -y/v\y-M + ^ 02 , (49) 



with 



fa = -5A 3 /(y - l) 4 + 3A(4r/ - A - A 2 )/(y - l) 3 + (4t?A - 4r? + A 

+7A 2 + A 3 )/(y - l) 2 + (-1 + 4??A - 28 V + 15A - A 2 - A 3 )/(y - 1) 

-1 + Ylyri - \\y\ + 7y - y 2 + 12r]\ - 24r/ + 14A - 11A 2 , (50) 

$ = ^J y 2 _ Ar} [ l5A 2 ^/(?/ - l) 3 + (10r?A£ 2 - 16^ + 24A£ - 10A£ 2 

-20A - 6A 2 0/(y - I) 2 + (-4 - 14r?A£ 2 - 48^ + 66r?£ 2 - 24r/£ 3 

+8r/ 2 £ 3 - 76A£ + 14A£ 2 + 48A + 3A 2 £ + 25£ - 26£ 2 + 8^)/(y - 1) 

+3 - 3y + 57r?£ - 22r?£ 2 - 12A£ - 21A + 34£ - 18£ 2 + 8£ 3 ] , (51) 

0f = 15(\ 2 /(y-l) 3 + (60r7CA-1677C- 3077C 2 A-16CA-21CA 2 

+ 10( 2 \)/(y - l) 2 + (-104?7CA - 84t]( + 52r]( 2 X + 122r]( 2 - 40r]( 3 
+160r/ 2 C - 160r/ 2 C 2 + 40r/ 2 C 3 - 24(A + 9(A 2 + 17C - 4( 2 A - 22( 2 
+8( 3 )/(y - 1) + 18 - 29yr] + 27y\ - 21y + 3y 2 + 46r/CA - 59r?C 
-14r^C 2 A + 78r?C 2 - 16r/C 3 - 71r/A + 59r? - AZrfC + 26r/ 2 C 2 - 8r/ 2 C 3 
+46r/ 2 - 6(A - 3CA 2 + 69( - 6( 2 A - 52( 2 + 16C 3 - 42A + 24A 2 , (52) 

where 

e = 2-C, X = p + V • (53) 
The function \1/ can be written in the form, 

_ f arccos to min - arccos u max , y < 1 , . 

[ arcosh u; max - arcosh u; min , y > 1 ' 

with 

_ 2(y - l)t mjn[maa; + y(y m -y)-2rj 

^min,max , • v""/ 

Due to terms containing inverse powers of (y — 1) the expression fl49| ) for Af(y) is 
apparently divergent at y — 1. However, expanding Af(y) in powers of (y — 1) for 
y < 1 and y > 1 one can check that this function is regular at y = 1. 

The HQET contribution to the decay distributions is known to render them un- 
reliable near the endpoint values of the tauon energy. This ambiguity reveals itself 
in the polarization as well. Similar problems appear also in calculations of perturba- 
tive corrections [16|, fl7l| . All these problems are cured if instead of distributions their 
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Figure 1: Integrated polarization of r lepton along the direction of virtual W in the 
Born approximation (dashed) and including HQET corrections (solid) as functions of 
the scaled r energy y. The mass of the b quark taken to be 4.75 GeV, c quark 1.35 
GeV. 

moments are considered |18|, |i~9| , |T0|| . In the case of r polarization a better defined 
quantity is the integrated polarization 

P (y)= T dy ( — - — ) I r dy ( — + — ) (56) 

where both the lowest-order perturbative and the HQET terms are included. In 
Fig.|T| the integrated polarization is shown as a function of the scaled energy y of the r 
lepton. The lowest order prediction corresponds to the dashed line and the solid line is 
obtained including HQET corrections. The question arises whether the perturbative 
QCD corrections can change this result significantly. As already suggested in the 
Introduction, it is plausible that no such thing happens. 

On integration over the whole range of the charged lepton energy one arrives at 
the total polarization at the tree level corrected for the 0(1 /ml) effects as predicted 
by HQET. For m& = 4.75 GeV and m c = 1.35 GeV, we obtain 

K G 

P = -0.7235 + 4.21— £ + 1.48— £ . (57) 
m{ mf 
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Taking K b = 0.15 GeV 2 and G b = -0.18 GeV 2 we obtain P = -0.706. 

Although we are mostly concerned with the tau lepton polarization here, the 
formulae derived in the present work may well be used in evaluating the polarization 
of the light leptons. Interestingly, in the limit of a vanishing mass of the charged 
lepton the polarization falls to zero apart from the endpoints, c.f. ( [49|) and (^). It 
is due to the chiral V — A structure of the weak charged current that, according to 
Eq. fllTD , the decay widths with a definite polarization differ by a term proportional 
to m T s M . The polarization four-vector of the charged lepton can be decomposed as 
follows: 

a" = (V, s) = (^/l - (s ± y, s ± , ^1 - ft.) 2 ) , (58) 

where s± is understood to mean the part of the three- vector s perpendicular to the 
direction of the charged lepton. The quantities E and p denote, respectively, the 
energy and the three-momentum value of the charged lepton. This form can easily 
be seen to meet the definition of the polarization four-vector, c.f. Appendix A. As 
the lepton mass approaches zero Eq.(^) gives 

m/ « sjl - (s ± ) 2 + m (0, s±, 0) . (59) 

However, if we want to keep the angle subtended by the polarization vector and the 
lepton momentum constant the parallel part of the polarization should be proportional 
to the perpendicular one, thereby forcing the factor of \Jl — (sj_) 2 to be of order 
of m/E. Then the r.h.s of Eq.(|5*9]) tends to zero for m — > 0. For the vanishing 
charged lepton mass the polarization can be non-zero only where the virtual W boson 
is collinear with the charged lepton. In general the contribution to polarization is 
appreciable only for W direction within the cone defined by the condition 

\s ± \/\s\=0(m/E) . (60) 

In particular this happens if p is much larger than the energy of the neutrino. For 
semitauonic B decays the condition (|60|) is satisfied in the whole phase space and the 
resulting polarization is fairly large. 

6 Summary 

The polarization of the tau lepton along the W boson direction in semileptonic B 
decays has been found at the tree level in perturbative QCD and the leading order 
HQET corrections have been included. The quantity is of experimental interest. The 
fact that it does but slowly vary in the regime of low energies of the charged lepton 
is rather favorable in this context ||. The QCD one- loop corrections are unknown 
but their irrelevance for the longitudinal polarization both in the rest frame of the 
decaying quark and that of the W boson indicates that no great change is to be 
expected once they are incorporated. 
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A The HQET calculations 

We will presently construct the four-vector s representing a charged lepton polarized 
along the direction of the W boson. The defining properties of s are 



and 



s 2 = -1 (61) 



s ■ r = , (62) 



complemented by the relation s || W. Since we are working in the rest frame of the 
decaying meson, the four-vector s can be decomposed as 

s = Av + BW, (63) 

where v and W stand for the four-velocity of the B meson and the four-momentum 
of the intermediate W boson, respectively. While this form automatically satisfies 
s || W, the other two relations ( j61|) and ( |62"D have to be imposed, hence yielding the 
expressions for the coefficients appearing in fl63l) . With v = (1,0,0,0), one readily 
identifies: 

v-r = y/2, v-u = x/2, v 2 = -1 . (64) 
These combined with the other dot products lead to the following formulae: 

A± = , (65) 

y/y{t + 7j)(y + x)-yH-(t + ri)* 

B± = . ±V ■ (66) 

yy(t + v){y + x ) - y 2t - (* + v) 2 

The evaluation of the HQET corrections involves differentiation over the neutrino 
energy, once or twice. The denominator in the above expressions is easily seen to 
vanish at the point where the W boson is at rest. It is known that the kinematics 
of the process, together with the delta functions and their derivatives, finally reduces 
to integration over the partonic phase space. Then there is one such point where the 
denominator vanishes, 

y = l-v/p + -^— , (67) 

i - VP 

t = (1-v/p) 2 • (68) 
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One might thus raise the question of analyticity of the expressions obtained in this 
way. However, the divergences cancel and moreover the resulting distribution is con- 
tinuous if we ignore the endpoint behaviour. That this is so indeed, may be verified 
by changing the variables from t to the square of the three-momentum of the W 
boson. Then the singularity makes its presence only on integration over w% rather 
than affecting the analytical structure of the distributions. It turns out that using 
this variable one obtains an analytic expression. This is made clear once one notices 
that the only terms that occur in the course of the calculation are the dot products 
of the four-vector s and the other four-vectors. Writing them out explicitly, 

s +- v = i 2 f fi ' ( 69 ) 

V2T - r 3 cos 2 

s+ . w= (y + f*™t>-^ , (70) 

l\/y — r 3 cos^ 

with 

cos^-^-^ 4 , (71) 

where the subscript denotes the polarization direction, we easily verify that the triple 
differential distribution is analytic in the neutrino energy. Lastly, let us note that 
another change of variable can be useful for evaluating the distribution. Namely, using 
the cosine of the angle subtended by the tau lepton and the neutrino jTT| eliminates the 



singular terms from the double differential distribution. We have checked numerically 
that the resulting distribution is the same. 
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